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Abstract 

We consider a multidimensional reaction-diffusion equation of either ignition or monos¬ 
table type, involving periodic heterogeneity, and analyze the dependence of the propagation 
phenomena on the direction. We prove that the (minimal) speed of the underlying pulsating 
fronts depends continuously on the direction of propagation, and so does its associated profile 
provided it is unique up to time shifts. We also prove that the spreading properties [24] are 
actually uniform with respect to the direction. 
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1 Introduction 

In this work, we focus on the heterogeneous reaction-diffusion equation 

d t u = div(A(x)Vit) + q(x) ■ Vu + f(x, u), M. N . (1) 

Here A = (Aij)i<ij<iy is a matrix field, and q = (qi, ..., gjv) is a vector field, to be precised later. 
The nonlinearity / is of either the monostable type (not necessarily with the KPP assumption) or 
ignition type, which we will define below. We would like to understand, in the periodic framework, 
how the propagation phenomena depend on the direction. 

On the one hand, we prove that the minimal (and, in the ignition case, unique) speed of 
the well known pulsating fronts depends continuously on the direction of propagation. On the 
other hand, we prove that the spreading properties are in some sense uniform with respect to the 
direction, thus improving the seminal result of Weinberger [24]. While in the ignition case, these 
properties will mostly follow from the well known uniqueness of the pulsating traveling wave, such 
uniqueness does not hold true in the monostable case where the set of admissible speeds is infinite. 
Our argument will be inspired by [10], [4], and will rely on an approximation of the monostable 
nonlinearity by some well-chosen ignition nonlinearities. 

1.1 Main assumptions 

Let Li,...,Ln be given positive constants. A function h : R N —> R is said to be periodic if 

h(x i,..., Xk + L k , ...,x N ) = h{xi, ...,x N ), 

for all 1 < k < N, all (x\ ,...,xn) € M. N . In such case, C = (0,Li) x ••• x (0 ,Ljv) is called the 
cell of periodicity. Through this work, we put ourselves in the spatially periodic framework and 
assume that 

for all 0 < i, j < N, the functions Ajj : —> R, : 1^ —> R are periodic, 

for all u £ R+, the function f(-,u) : R w —> R is periodic. 

Moreover, we assume that A = {Aij)i<i j<jy is a C 3 matrix field which satisfies 

A(x) is a symmetric matrix for any x £ R^, 

30 < cn < a 2 < oo, V(x,0 £ R N x R N , <xi|^| 2 < J2i,j < a 2 |C| 2 - 

Concerning the advection term, we assume that q = (q±, ..., qat) is a C 1,s vector field, for some 
6 > 0, which satisfies 

div <7 = 0 in R w and VO < i < N, J = 0. (4) 

The advection term in the equation is mostly motivated by combustion models where the dynamics 
of the medium also plays an essential role. In such a context, the fact that the flow q has zero 
divergence carries the physical meaning that the medium is incompressible. 

Furthermore, we will assume that / satisfies either of the following two assumptions. 

Assumption 1.1 (Monostable nonlinearity). The function f : R^ x R + — y R is of class C 1,a in 
( x,u ) and C 2 in u, and nonnegative on R w x [0,1]. Concerning the steady states of the periodic 
equation (1), we assume that 

(i) the constants 0 and 1 are steady states (that is, /(•, 0) = /(•, 1) = 0 in R N ); 

(ii) Vu £ (0,1), 3x £ R^, f(x, u) > 0. 

(Hi) there exists some p > 0 such that f(x,u) is nonincreasing with respect to u in the set 
R N x (1-/»,!]. 
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Notice that, if 0 < p(x) < 1 is a periodic stationary state, then p = 0 or p = 1. Indeed, since 
f{x,p) > 0, the strong maximum principle enforces p to be identically equal to its minimum, thus 
constant and, by (ii), the constant has to be 0 or 1. 

Assumption 1.2 (Ignition nonlinearity). The function / : K w x R + —> R is locally Lipschitz- 
continuous on x R + . Concerning the steady states of the periodic equation (1), we assume 
that 

(i) there exists 0 < 9 < 1 such that 

VO < u < 9, Vs e M. N , /( x, u) = 0, 

as well as 

Vx&R N , f(x, 1) = 0. 

(ii) Vit € (0,1), 3a; £ R N , /( x, u) > 0. 

(in) there exists some p > 0 such that f(x,u) is nonincreasing with respect to u in the set 
R N x (1 — p, 1]. 

Notice that, similarly as above, (ii) implies that if 6 < p(x) < 1 is a periodic stationary state 
then p = 6 or p = 1 . 

1.2 Comments and related results 

Under Assumption 1.1, Assumption 1.2, equation (1) is referred to as the monostable equation, 
the ignition equation respectively. Both sets of assumptions arise in various fields of physics and 
the life sciences, and especially in combustion and population dynamics models where propaga¬ 
tion phenomena are involved. Indeed, a particular feature of these equations is the formation of 
traveling fronts, that is particular solutions describing the transition at a constant speed from one 
stationary solution to another. Such solutions have proved in numerous situations their utility in 
describing the spatio-temporal dynamics of a population, or the propagation of a flame modelled 
by a reaction-diffusion equation. 

Equation (1) is a heterogeneous version of the well known reaction-diffusion equation 

d t u = Au + f(u), (5) 

where typically / belongs to one of the three following classes: monostable, ignition and bistable. 
Homogeneous reaction-diffusion equations have been extensively studied in the literature (see [16], 
[2, 3], [11], [9], [23] among others) and are known to support the existence of monotone traveling 
fronts. In particular, for monostable nonlinearities, there exists a critical speed c* such that all 
speeds c > c* are admissible, while in the bistable and ignition cases, the admissible speed c = c* 
is unique. Moreover, in both cases, the speed c* corresponds to the so-called spreading speed of 
propagation of compactly supported initial data. 

Among monostable nonlinearities, one can distinguish the ones satisfying the Fisher-KPP as¬ 
sumption, namely u K > is maximal at 0 (meaning that the growth per capita is maximal 

at small densities), the most famous example being introduced by Fisher [12] and Kolmogorov, 
Petrovsky and Piskunov [17] to model the spreading of advantageous genetic features in a popu¬ 
lation: 

dtU = A u + u(l — u). 

Let us notice that our work stands in the larger class of monostablc nonlinearities. 

Nevertheless, much attention was more recently devoted to the introduction of some hetero¬ 
geneity, taking various forms such as advection, spatially dependent diffusion or reaction term. 
Taking such a matter into account is essential as far as models are concerned, the environment 
being rarely homogeneous and may depend in a non trivial way on the position in space (patches, 
periodic media, or more general heterogeneity...). We refer to the seminal book of Shigesada and 
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Kawasaki [22], and the enlightening introduction in [7] where the reader can find very precise and 
various references. As far as combustion models are concerned, one can consult [9], [25] and the 
references therein. 

Traveling front solutions in heterogeneous versions of (5) with periodicity in space, in time, or 
more general media are studied in [24], [10], [15], [26], [4], [8], [19], [20] among others. For very 
general reaction-diffusion equations, we refer to [5] for a definition of generalized transition waves 
and their properties. 

In this work, we restrict ourselves to the spatially periodic case, which provides insightful 
information on the role and influence of the heterogeneity on the propagation, as well as a slightly 
more common mathematical framework. In this periodic setting, let us mention the following 
keystone results for ignition and monostable nonlinearities. Weinberger [24] exhibited a direction 
dependent spreading speed for planar-shaped initial data and proved, in the monostable case, 
that this spreading speed is also the minimal speed of pulsating traveling waves moving in the 
same direction. His approach relies on a discrete formalism, in contrast with the construction of 
both monostable and ignition pulsating traveling waves by Berestycki and Hamel [4], via more 
flexible PDE technics. In this PDE framework, note also the work of Berestycki, Hamel and 
Roques [8] where KPP pulsating fronts are constructed without assuming the nonnegativity of the 
nonlinearity. Our main goal is to study how these results behave when we vary the direction of 
propagation. 

Let us give another motivation for our analysis of the dependence of the propagation on the 
direction. Our primary interest was actually to study the sharp interface limit e —0 of 

9 t u e = eAu £ + i/(J,u £ ), (6) 

arising from the hyperbolic space-time rescaling u e (t,x) := u (|, |) of (1), with A = Id, q = 0. 
The parameter e > 0 measures the thickness of the diffuse interfacial layer. As this thickness 
tends to zero, (6) converges — in some sense— to a limit interface, whose motion is governed 
by the minimal speed (in each direction) of the underlying pulsating fronts. This dependence of 
the speed on the (moving) normal direction is in contrast with the homogeneous case and makes 
the analysis quite involved. In particular, it turns out that we need to improve (by studying the 
uniformity with respect to the direction) the known spreading properties [24], [4], for both ignition 
and monostable nonlinearities in periodic media. We refer to [1] for this singular limit analysis, 
using some of the results of the present work. 

2 Main results 

Before stating our main results in subsection 2.2, let us recall the classical results on both pulsating 
fronts and spreading properties in subsection 2.1. 

2.1 Pulsating fronts and spreading properties: known results 

The definition of the so-called pulsating traveling wave was introduced by Xin [25] in the framework 
of flame propagation. It is the natural extension, in the periodic framework, of classical traveling 
waves. Due to the interest of taking into account the role of the heterogeneity of the medium on the 
propagation of solutions, a lot of attention was later drawn on this subject. As far as monostable 
and ignition pulsating fronts are concerned, we refer to the seminal works of Weinberger [24], 
Berestycki and Hamel [4]. Let us also mention [8], [13], [14], [19] for related results. 

For the sake of completeness, let us first recall the definition of a pulsating traveling wave for 
the equation (1), as stated in [4]. 
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Definition 2.1 (Pulsating traveling wave). A pulsating traveling wave solution, with speed c > 0 

is an entire solution u(t,x ) — t £ M, x £ M . N — of (1) satisfying 


in the direction n £ 1 


N 


Vfc £ LjZ, 


i k ■ n 
u(t, x) = u I t H-, x + k 


for any t £ R and x £ M. N , along with the asymptotics 

u(— oo, •) = 0 < u(-, •) < u(+oo, •) = 1, 

where the convergences in ±oo are understood to hold locally uniformly in the space variable. 

One can easily check that, for any c > 0 and n £ S N ~ X , u(t,x) is a pulsating traveling wave 
with speed c in the direction n if and only if it can be written in the form 

u(t, x) = U(x ■ n — ct, x), 

where U(z,x) — z £ R, x £ R N — satisfies 


for all z £ R, U{z, ■) : 


i,N 


is periodic, 


U{— oo, •) = 1 <[/(•,•) < U(+oo, •) = 0 uniformly w.r.t. the space variable, 
along with the following equation 

di v x (AX7 x U) + (n • An) d zz U + di v x (And z U) + d z (n ■ AX7 X U) 

+ q ■ V X U + (q ■ n) d z U + cd z ll + /( x, U) = 0, on M x R N . 
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We can now recall the results of [4] (see also [24] for the monostable case), on the existence of 
pulsating traveling waves for the spatially periodic monostablc and ignition equations. Precisely, 
the following holds. 

Theorem 2.2 (Monostable and ignition pulsating fronts, [4],[24]). • Assume that f is of the 

spatially periodic monostable type, i.e. f satisfies (2) and Assumption 1.1. Then for any 
n £ there exists c*{n ) > 0 such that pulsating traveling waves with speed c in the 

direction n exist if and only if c> c* (n). 

• Assume that f is of the spatially periodic ignition type, i.e. f satisfies (2) and Assump¬ 
tion 1.2. Then for any n £ S w_1 , there exists a unique (up to time shift) pulsating traveling 
wave, whose speed we denote by c*(n ) > 0. 

Furthermore, in both cases, any pulsating traveling wave is increasing in time. 

The introduction of these pulsating traveling waves was motivated by their expected role in 
describing the large time behavior of solutions of (1) for a large class of initial data. In this 
context, let us state the result of [24] for planar-shaped initial data. 

Theorem 2.3 (Spreading properties, [24]). Assume that f is of the spatially periodic monostable 
or ignition type, i.e. f satisfies (2) and either of the two Assumptions 1.1 and 1.2. Let uq be a 
nonnegative and bounded initial datum such that ||ito ||oo < 1 and 

3C >0, x ■ n > C => uo(x) = 0, 

liminf uq(x) > 0 (monostable case), liminf uq(x) > 6 (ignition case) 


x-n—t—oo 
5JV-1 


X'Tl—y — OO 


for some n £ Er 

Then the solution u of (1), with initial datum uq, spreads with speed c*(n) in the n-direction 
in the sense that 

Vc < c*(n), lim sup |1 — u(t, x)\ = 0, (8) 


Vc > c*(n), 


lim sup |1 — u(t, x)\ = 0, 

t->+oo x . n < c t 

lim sup u(t,x) = 0. 

t-H-oo x . n>ct 


(9) 



Remark 2.4 (Link between spreading speed and wave speed). In [24], Weinberger was actually 
concerned with a more general discrete formalism where pulsating waves are not always known to 
exist. Therefore, the fact that the spreading speed and the minimal wave speed are one and the 
same was only explicitly stated in the monostable case. 

However, under the ignition Assumption 1.2 and benefiting from the results in [4], it is clear 
by a simple comparison argument that the solution associated with any such initial datum spreads 
at most with the wave speed c*(n), namely (9) holds true. Furthermore, one may check, using for 
instance U*(x ■ n — (c*(n) — a)t,x ) ~ 5 as a subsolution of (1), where U* is the pulsating wave 
with speed c*(n) and a > 0, 5 > 0 are small enough, that (8) also holds true, at least for some 
large enough initial data. Thus, the spreading speed exhibited by Weinberger must be c*{n), as one 
would expect. 

We will use a very similar argument in Section 1, which is why we omit the details. Moreover, 
it is a simplification of a classical argument, which originates from [11] in the homogeneous frame¬ 
work, and usually aims at proving the stronger property that the profile of such a solution u(t , x) 
of the Cauchy problem converges to that of the ignition pulsating wave. We refer for instance 
to the work of Zlatos [27], which dealt with a fairly general multidimensional heterogeneous (not 
necessarily periodic) framework, and covers the above result under the additional assumption that 
f(x,u ) is bounded from below by a standard homogeneous ignition nonlinearity. 

Various features of pulsating fronts and many generalizations of spreading properties have been 
studied recently. Nevertheless, as far as we know, nothing is known on the dependence of these 
results on the direction of propagation. Our results stand in this new framework and are stated 
in the next subsection. 

2.2 Pulsating fronts and spreading properties: varying the direction 

As recalled above, the periodic ignition equation admits a unique pulsating traveling wave in any 
direction n £ while the periodic monostable equation (1) admits pulsating traveling waves 

in any direction n £ S^ -1 , for any speed larger than some critical c*(n ) > 0. The latter is a 
consequence of the former, as was proved in [4] by approximating the monostable equation with 
an ignition type equation. With some modifications of their argument, we will prove the following 
continuity property. 

Theorem 2.5 (Continuity of minimal speeds). Assume that f is of the spatially periodic monos¬ 
table or ignition type, i.e. f satisfies (2) and either of the two Assumptions 1.1 or 1.2. 

Then the mapping n £ S N_1 H > c*(n) is continuous. 

In the Fisher-KPP case the continuity of the velocity map n c*(n), even if not explicitly 
stated, seems to follow from the characterization of c*(n) (see [24], [4]). However, for other types 
of nonlinearities (and in particular, in the more general monostable case), such a property seems 
to be far from obvious. 

For the sake of completeness, let us state the continuity of the profile of the ignition wave, 
which will be proved simultaneously. 

Theorem 2.6 (Continuity of ignition waves). If f satisfies (2) and the ignition Assumption 1.2, 
then the mapping 

n £ i->- U*(z,x;n) 

is continuous with respect to the uniform topology, where 

u*(t,x;n ) = U* (x ■ n — c* ( n)t , x; n) 

is the unique pulsating traveling wave in the n direction, normalized by min,j, eR iv U* (0, x; n) = 

In Section 3, we deal with the ignition case, proving both the continuity of the speed (Theorem 
2.5) and that of the profile (Theorem 2.6). To do so we take advantage of the uniqueness of the 
pulsating wave in each direction. 
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Then, in Section 4, we approach our original monostable equation by some ignition type prob¬ 
lems, and prove that the associated ignition speeds converge to c*(n) not only pointwise (as in [4]), 
but even uniformly with respect to n € S^ -1 . The continuity of the minimal speed (Theorem 2.5) 
then immediately follows. Unfortunately, the lack of a rigorous uniqueness result of the monos¬ 
table pulsating wave with minimal speed (at least up to our knowledge) prevents us from stating 
continuity of its profile with respect to the speed of propagation. We refer to [14] for uniqueness 
results in the Fisher-KPP case and discussion on the general monostable framework. 

We also stated above the well known fact that for any planar-like initial data in some direc¬ 
tion n, the associated solution of (1) spreads in the n direction with speed c*(n). Our main result 
consists in improving (compare Theorem 2.7 with Theorem 2.3) this property by adding some 
uniformity with respect tons as follows. 


Theorem 2.7 (Uniform spreading). Assume that f is of the spatially periodic monostable or 
ignition type, i.e. f satisfies (2) and either Assumption 1.1 or Assumption 1.2. Let a family of 
nonnegative initial data (uo, n )raeS JV - 1 be such that 


3C > 0, Vn6§ w \ x ■ n > C => uo, n (x) = 0, 

3 p > 0 (ignition case) 

3p > 0 ( monostable case ) 

sup sup Uo,n(%) < 1. 
nes^- 1 xeu" 

We denote by (u n ) ne §N-i the associated family of solutions of ( 1 ). 

Then, for any a > 0 and 6 > 0, there exists r > 0 such that for all t > t, 


3K > 0, inf uq n (x ) > u, 

nGS-W -1 , x-n<—K 


( 10 ) 

(11) 

( 12 ) 


sup sup |1 — u n (t, x)\ < 6 , (13) 

nGS ^ -1 x-n<(c* (n)—a)t 

sup sup u n (t,x) < <5. (14) 

nG S N_1 x-n>{c*{n)-\-Oi)t 

The difficult part is again to deal with the monostable case. The proof of the lower spreading 
property (13) will again rely on an ignition approximation of the monostable equation, whose 
traveling waves will serve as nontrivial subsolutions of (1). This is performed in Section 5. Then, 
Section 6 is devoted to the proof of the upper spreading property: we prove (14) in subsection 6.1 
and, for sake of completeness, relax assumption (12) in subsection 6.2. 

Last, in Section 7, we prove Theorem 2.7 in the ignition case. 


3 Continuity of ignition waves 

Let us here consider a periodic nonlinearity / of the ignition type, namely satisfying Assump¬ 
tion 1.2. As announced, we will prove simultaneously the continuity of both mappings n H > c*(n) 
and n H > U*(z,x\n ), where we recall that c*{n) and U*(x ■ n — c*{n)t,x\n) denote respectively 
the unique admissible speed and the unique pulsating wave in the direction n , normalized by 

min U*(0,x-n) = (15) 

Proofs of Theorem 2.5 (ignition case) and Theorem 2.6. We first claim (we postpone the proof to 
the end of this section) that 

k \= inf c*(n) > 0. (16) 

neS "- 1 

Let us now prove that that c*(n) is also bounded from above, using 
(t, x) ^ v(t, x) := min{l, 0 + Ce~ x{x - n - 2aiXt) } 
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as a supersolution of (1). Here C and A are positive constants to be chosen later, and aq comes 
from hypothesis (3). Indeed, when v < 1, it satisfies 

dtv — div (A(x)Vzi) — q(x ) • Vu — f(x, v) 

= [2aiA 2 - (n ■ An) A 2 + A div (An) + A q-n) x Ce~ x ^ n - 2aiXt) - /(x, v) 

> [aiA 2 - A |div (An )| - A|g • n\ - M] x Ce~ x ( xn ~ 2aiXt) > 0, (17) 

where 

M := sup ^^<+oo (18) 

comes from the Lipschitz continuity of /, and the last inequality holds provided that A is large 
enough, independently of n £ S N ~ 1 . As 1 is a solution of (1), it is then clear that v is a generalized 
supersolution of (1). Then, choosing C > 0 so that v(t = 0, x) lies above the traveling wave 
u*(t = 0, x; n) = U*(x ■ n, x ; n) at time 0, we can apply the comparison principle and obtain that 
c*(n) < 2a±\. Putting this fact together with (16), we conclude that 

0 < k := inf c*(n) < sup c*(n) =: K < +oo. (19) 

neS"- 1 neS'V” 1 


We now let some sequence of directions nk —> n S § w_1 . As we have just shown, the sequence 
c*(nk) is bounded and, up to extraction of a subsequence, c*(nk) —> c > 0. We also choose 
the shifts Zk so that, for all k, max iel » U*(zk, x; nk) = 0. In particular, recalling that U* is 
monotonically decreasing with respect to its first variable, we have for all k that 

Vz > £fc,Vx £ K^, 0 < U*(z,x;nk) < 0. 


Then 

u k (t,x) := U*(z k + x ■ n k ,x + c*(n k )tn k ] n k ) 

satisfies 

d t u k = div (A(x)Vu k ) + q(x) ■ Vu k + c*(n k )Vu k ■ n k , (20) 

for all t £ R and all x in the half-space x-n k > 0 (recall that U* solves (7) and that, in the ignition 
case, /(x, u) = 0 if 0 < u < 9). 

Let us now find a supersolution of (20) of the exponential type, namely 

u k (t,x) := (j> k (t,x) x e~ x ° x ' nk , (21) 


where cf> k will be a well-chosen positive and bounded function. 

For any n £ one may define (see Proposition 5.7 in [4]) the principal eigenvalue problem 

J -L n} x4> n ,\ = n(n, X)(/)n,x in R^, 

) 4>n,x > 0 is periodic, 


where 


L n ,\4> div (AVtfi) + A 2 (n ■ An)<f> — A(di v(An(f>) + n ■ AV0) + q ■ Vtp — A (q ■ n + k)</>, 
with k > 0 given by (16). In the sequel, the eigenfunction <p n ,\ is normalized so that 

min (f>n t \(x) = 9. 

xEC 

As stated in Proposition 5.7 of [4], the function A >->■ /i(n, A) is concave and satisfies, for any n, that 
fi(n, 0) = 0 (any positive constant is clearly a principal eigenfunction of —L nj o), and d\fj,(n,0) = 
k > 0 . 



It follows that one can find some small Ao > 0 such that, for any n £ S w 1 , 

H(n, A 0 ) > 0. 

Indeed, proceed by contradiction and assume that for any j £ N*, there exists rij such that 
fi(nj,l/j) < 0. Then, by y(rij, 0) = 0 and by concavity, one has that y(rij, A) < 0 for all A > 4. 
By uniqueness of the principal normalized eigenfunction, it is straightforward to check that y(n, A) 
depends continuously on both n and A, as well as with respect to the uniform topology. Thus, 
one can pass to the limit and conclude that y(noo, A) < 0 for some noo = lim rij (up to extraction 
of a subsequence) and all A > 0. This contradicts the fact that d\/i(noo, 0) = n > 0. 

Notice that, by continuity of the eigenfunction with respect to n and A in the uniform topology, 
it is clear that for any bounded set A, 

max max max<A„ \(x) < +oo. (22) 

neS^-i AeA 16 c ’ 


Choosing Ao as above and 


<t> k {t,x) := 4>n k ,\ 0 {x + c*(n k )tn k ), 

in (21), one gets that 

d t u k - div (A(x)Vu/c) ~ q{x) ■ Vu k - c*(n k )Vu k ■ n k 
= [c*(n k )n k • V</) nfci A 0 - L nkt \ 0 cf> k + A 0 (-k + c*{n k ))<j> k - c*(n k )n k ■ V</> nfci A 0 ] x e _Ao(x ' nfc) 
= [p(n k , A 0 ) + A 0 (c*(n fc ) - n)\u k > 0. 

In other words, as announced, u k is a supersolution of (20). 

Let us now prove that 


V£ el, Vx ■ n k > 0, u k (t,x) < u k {t,x). (23) 

Proceed by contradiction and define a sequence of points such that 

u k (tj , Xj ) u k (tj , Xj ) sup (u k (t,x) - u k {t,x)) > 0. 

tGl,x-nfe> 0 

Now write Xj = (xj ■ n k )n k + yj for any j > 0. Note that, since u k (t,x ) and u k (t,x) both tend 
to 0 as x ■ n k —> +oo uniformly with respect to t, then (xj ■ n k )j^N must be bounded. Thus, 
up to extraction of a subsequence, we can assume that Xj ■ n k —> a > 0 as j —> oo. Moreover, 
since yj is orthogonal to n k , since <f>n k ,\ 0 is periodic and since U* is periodic with respect to its 
second variable, we can assume without loss of generality that yj + c*{n k )tjn k € C the cell of 
periodicity. As yj is orthogonal to n k for all j £ N, we can extract a subsequence such that both 
tjj — y y 0 o £ and tj — y £ M. 

Finally, u k — u k reaches its positive maximum, over t £ R and x ■ n k > 0, at (t = too, x = 
an k + yoo)- Moreover, as 

Vz • n k = 0, i/fc(0, x) <6 < u k (0, x), 

the maximum is reached at an interior point, which contradicts the parabolic maximum principle. 
Thus, (23) is proved. 

Now, by standard parabolic estimates and up to extraction of a subsequence, we can assume 
that, as k —y oo, the sequence u* (t — , x; n^j = U*(x ■ n k — c*{n k )t + z k ,x;n k ) converges 

locally uniformly, along with its derivatives, to a solution Uoo(t,x ) of (1). Moreover, u 0 0 satisfies 

VZ £ n^LiZ, Uoo{t,x) = Uoo (t + —— ,x + l 
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In a similar way than the discussion after Definition 2.1 of pulsating waves, this means that 
Uoo(t,x) = Uoo(x ■ n — ct,x) where U 00 (z,x) is periodic with respect to its second variable and 
satisfies 

dxv x (A\7 X U) + (n ■ An) d zz U + div x (And z U) + d z {n ■ AV X U) 

+ q ■ S7 X U + (q ■ n) d z U + cd z U + f(x, U) = 0 on R x R w . 

It is then straightforward to retrieve that the sequence U* (z + z k , x; nk) also converges, along with 
its derivatives, to this function U 00 (z,x). In particular, Uoo is nonincreasing with respect to its 
first variable, and satisfies the inequalities 

0 < Uoo{z,x) < 1 . 

Furthermore, noticing that u* (t — e *(n fc ) , x\ n^j = u k (t, x — c*(n k )tn k ), it follows from passing to 
the limit in (23), and thanks to (22), that 

u 00 (t,x)<Ae- x °( x - n ~ ct \ 


for some A > 0 and all x ■ n> ct. 

Thus, U oo (x ■ n — ct,x) < Ae~ x A x ' n ~ ct \ for all t £ R. and x ■ n> ct. This means that Uoo(z, x) 
converges exponentially to 0 as z —> +oo, uniformly with respect to its second variable: 

Vz > 0, Vx £ R w , Uoo{z,x) < Ae~ x ° z . (24) 

By monotonicity with respect to its first variable, Uoo(z, x) converges as z —> —oo to some periodic 
function p(x). Or, equivalently, Uoo{t,x) converges as t —> +oo to the same function p(x). By 
standard parabolic estimates, we get that p(x) is a periodic and stationary solution of (1). Let us 
show that p = 1. From our choice of the shifts Zk and up to extraction of another subsequence, 
there exists some Xoo such that Uoofi^Xoo) = 9 , hence maxp > 9. Assume first that maxp = 9. 
Then Uco(t, x) < 9 for all t€R and x £ R N and, by the strong maximum principle, Uoq = 9. This 
contradicts the inequality (24) above. Therefore, maxp > 9. Using again the strong maximum 
principle and the fact that f(x,u) = 0 for all u <9 and x £ M. N , we reach another contradiction 
if min p < 9. Therefore min p > 9 and, thanks to part (ii) of our ignition Assumption 1.2, p = 1 
the unique periodic stationary solution of (1) above 9. 

From the above analysis, we conclude that [/oo(-,-) = U*(- + Z,-',n) the unique pulsating 
traveling wave in the n direction with speed c = c*(n), where Z is the unique shift such that 
maXj-gRV U*{Z , x: n ) = 9. This in fact proves, by uniqueness of the limit, that the whole sequence 
c*(rife) converges to c*(n), and that the whole sequence U*(- + Zk, ■; Ufe) converges locally uniformly 
to U*(-+Z, •; n). This in particular shows the continuity of the map n i-» c*(n), that is Theorem 2.5 
in the ignition case. 

Let us now conclude the proof of Theorem 2.6. Let us first prove that the sequence of shifts 
Zk is bounded. The normalization (15) implies that U*(0,yk',nk) = -4p, for some yk £ C that 
(up to some subsequence) converges to some y £ C. Since U*(zk,yk',n-k ) —> U*{Z,y\n) < 9 and 
E/*(0, y,n) = ^ 2 ^, the monotonicity of traveling waves enforces Zk > 0 for k large enough. Now 
proceed by contradiction and assume that (up to some subsequence) Zk —> +oo. Then, for all 
-z k < z < 0, 

U*(z + z k ,y k \n k ) < U*(0,y k -,n k ) = 

Passing to the limit as k +oo, we get that 

1 _|_ Q 

C^oo(z,2/) < —, 

for all z < 0. This contradicts the fact that Uoo is a pulsating traveling wave and converges to 1 
as z —> — oo. 

From the boundedness of the sequence z k , we can now rewrite the convergence as follows: the 
sequence U*(-,-;n k ) converges locally uniformly to !/*(•,•;n). It now remains to prove that this 
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convergence is in fact uniform with respect to both variables. Note first that uniformity with 
respect to the second variable immediately follows from the periodicity. Furthermore, for a given 
S > 0, let K > 0 be such that, for any x £ W N , 

r c 

0 < U*(z,x\n ) < - and 1 — - < U*(—z,x\ri) < 1, for all z > K. (25) 

From the locally uniform convergence with respect to the first variable, we have, for any k large 
enough, 

\\U*(; •; n k ) — U*(-, •; n) ||l~([_a-,a-] x r) < ^ ■ 

In particular, U*(K, x; n k ) < S and 1 — <5 < U*{—K , x; n k )- so that, by monotonicity with respect 
to the first variable, for any x £ R^ and k large enough, 

0 < U*(z,x;n k ) < S and 1 — S < U*{—z,x\n k ) < 1, for all z> K. (26) 

Combining (25) and (26), we get 

II U*(-,-;rik) - t^*(-, S ^)l|L“((-oo,-rc)u(if,oo)xR) < 

for any k large enough. As a result the convergence of U*(-,-;n k ) to {/*(•,-;n) is uniform in 
R x R w . This ends the proof of the continuity of ignition waves, that is Theorem 2.6. □ 

Proof of claim (16). Proceed by contradiction and assume that there exists a sequence rik £ 
such that c*(rik) —> 0. 

Now for any k , recall that the pulsating wave is normalized by 

mil Sv U*(0, x\ n k ) = —. (27) 

igl*' 2 

Up to extraction of a subsequence, we can assume as above that n k —>• n and 

u* ( t,x;n k ) -»• Uoo(t,x), 

where the convergence is understood to hold locally uniformly, and Uoo(t,x) is a solution of (20). 
By the strong maximum principle, it is clear that 0 < u 0 0 < 1. We also know, by the monotonicity 
of [/*(■, -;n k ) with respect to its first variable, by (27) and by passing to the limit, that 

1 _|_ Q 

Uqo{ t, x) > —-—, Vx • n < 0. 

Applying Weinberger’s result (see Theorem 2.3 as well as Remark 2.4), we get that the solution 
spreads at least at speed c*(n). In particular, as t —> +oo, Uoo(t,x) converges locally uniformly 
to 1. 

On the other hand, we fix x £ R w and s > 0, then we let some vector l £ be such 

that l ■ n > 0. In particular, for any large k, one also has that l ■ n k > ^ > 0. Then, for all 
large k, using the fact that c*{n k ) —> 0 and the monotonicity of «*(•, •; n k ) with respect to its first 
variable, we have that 


u* (s, x; n k ) < u* 


l ■ n k 
c*(n k ) 


, x; n k 


By passing to the limit as k —> +oo, we obtain that 


u* (0, x — l; n k ). 


u oo(s,x) < 1X00(0, x-l) < 1 , 

for all x £ R w and s > 0. This contradicts the locally uniform convergence of Uoo(t,x) to 1 as 
t —> + 00 . The claim is proved. □ 
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4 Continuity of the monostable minimal speed 

Let us here consider a periodic nonlinearity / of the monostable type, namely satisfying Assump¬ 
tion 1.1. We will prove the continuity of the mapping n H > c*(n), that is Theorem 2.5. To do so, 
we introduce a family f e (x,u), for small £ > 0, of ignition nonlinearities which serve as approxi¬ 
mations from below of the monostable nonlinearity f(x,u). Our aim is to prove that, by passing 
to the limit as £ —> 0, we indeed retrieve the dynamics of the monostable equation. This will be 
enough to prove Theorem 2.5. 

The family (/ e ) e , for small enough £ > 0, is chosen as follows: 

( v« e [-e, o], f e {x,u) = o 

Va; e R^, < VuG[0,1-e], f e (x,u) = f(x,u) 

{ Vu G [1 - £, 1 - §] , f s (x, u) = f (x, 1 - £ + 2{U - (1 - £))) . 

Notice that || f e — /||l=°(- e ,i) —> 0 as e —► 0, and that, thanks to Assumption 1.1 (in), f e lies below 
/ and 0 < £ < e' implies f e > f e >. Also, the equation 

d t u = div(A(x)Vu) + q(x) ■ Vu + f e (x, u), (28) 

where u is to take values between — e and 1 — |, is of the ignition type in the sense of Assumption 1.2 
(where 0, 6, 1 are replaced by — e, 0 and 1 — | respectively). In particular, for each n G 
there exists a unique ignition pulsating traveling wave 

u*(t,x\n) = U*(x ■ n — c* (n)t, x; n) 

of (28) in the n direction with speed c*(n) > 0, normalized by 

min U*(0,x:n) = 

iei» v ’2 

Furthermore, we have already proved in the previous section that the mappings n i —> c*(n) and 
n U*(-, •; n ) are continuous (with respect to the uniform topology). 

Theorem 4.1 (Convergence of speeds). Assume that f is of the spatially periodic monostable 
type, i.e. f satisfies (2) and Assumption 1.1. Let f e (x,u) be defined as above. 

Then, as e ^ 0, c*(n) /* c*(n) uniformly with respect to n G 

As mentioned before, pointwise convergence was shown in [4], where the goal was to prove 
existence of monostable traveling waves for the range of speeds [c*(n), +oo). Here we prove that 
the convergence is actually uniform, which together with the continuity of speeds in the ignition 
case, immediately insures the continuity of n <—> c*(n), that is Theorem 2.5 in the monostable case. 

Proof. First note that, for any fixed n G and £ > 0, c*(n) < c*(n). Indeed, recalling that 

Uf (z, x; n) connects 1 — | to — e, one can find some shift ZGl such that XJ*(z + Z,x;n) < 
U*(z,x;n), where U* denotes a monostable pulsating traveling wave — connecting 1 to 0— with 
the minimal speed c*(n). By a comparison argument, it follows that c*(n) < c*(n). It is also very 
similar to check that, for any n G S^ -1 , 0 < £ < e' implies c*(n) > c*,(n). 

Let us now consider some sequences £*, —> 0 and nk —> n. Consider the estimate (19) where 
k and K should a priori depend on e. First, it is clear from the above that n(e) := inf n c*(n) is 
nonincreasing with respect to e. Also, since 

J ^ X , u\ 

sup M e := sup sup 1 — < +oo 

0<e<eo 0<e<eo x6R Ar , u£[—e,l —§] M 

(compare with (18)), arguing as we did to derive (19), we see that K(e) := sup n c*(n ) is uniformly 
bounded from above. As a result, we have 

0 < At := inf inf c*(n ) < sup sup c*(n) =: K < +oo. (29) 

0<e<E O neS N - 1 0<e<£o neS^- 1 
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Hence, we can assume, up to extraction of a subsequence, that c* k (n k ) —> c 00 > 0 as k —> oo. In 
order to prove Theorem 4.1, we have to prove that c 00 = c*{n). 

We begin by showing that U* (z, x; n k ) converges as k —> oo to a monostable pulsating traveling 
wave of (1), up to extraction of a subsequence. Indeed, proceeding as before, one can use standard 
parabolic estimates to extract a converging subsequence of pulsating ignition traveling waves, such 
that 

U* k (z,x;n k ) -x U oc (z, x), 

as k —> +oo locally uniformly with respect to (z,x) elx R w . Furthermore, 0 < U 00 (z,x) < 1 
solves (7) with c = Coo, is nonincreasing with respect to z, periodic with respect to x, and satisfies 
m i n a;eR N Uoo(0, x) = In particular, Uoo converges as 2 —> ±oo to two periodic stationary 
solutions of (1), which under the monostable Assumption 1.1 can only be 0 and 1. We can 
conclude that is a monostable pulsating traveling wave with speed Coo, hence Coo > c*(n). 

We now prove that Cqo = c*(n). Notice that f e lies below / but, since the direction varies, we 
cannot use a simple comparison argument to conclude that Coo < c*(n). Instead, we will use a 
sliding method as in [4]. To do so, we shall need the following lemma. 

Lemma 4.2 (Some uniform estimates). There exists C > 0 such that, for any small e > 0 and 
n £ the ignition pulsating traveling wave U* (z, x;n) satisfies 

\d zz U*(-, •; n)| < Cd z Uf (*, *, n), \V x d z U* e {-, ■; n)| < -Cd z U* e {;-,n). 

Proof. Let us define u*[t,x) := U*{x ■ n — c*(n)t,x;n). Then v(t,x) := d t u*(t,x) > 0 satisfies 

d t v = div(A(a;)Vu) + q(x) ■ Vu + v d u f e (x, u*), a.e. in R x l‘ ¥ . 

From our definition of the ignition approximation f e (x,u), we see that ||II§)) is 
uniformly bounded, independently on small e > 0 and n £ Therefore, from the interior 

parabolic L p -estimates (see [21, Theorem 48.1] for instance) and Sobolev embedding theorem, one 
gets 

\/(t 0 ,x 0 ) elx R n , |9 t u(to,xo)| + |V x w(< 0 , £o)| < Ci max v(t,x), (30) 

to — »I 3 ?—^ol^l 

for some C\ > 0 which is independent on to, xq , small e > 0 and n £ E> N ~ 1 . 

Furthermore, for any n £ S^ 1 , choose N integers fcj(n) £ {—1,0,1} such that 

k(n)L-n= max (k\L\, ...IznLn) • n, 

where k{n)L := [k\{n)Li, ...kif{n)Lifi). Then 

0 < inf k{n)L ■ n < sup k{n)L ■ n < +oo, 
neS^- 1 ne S^- 1 


and hence, thanks to (29), 


0 < inf 

0 <£<£o,nG § Ar_1 


k(ri)L • n 

c e( n ) 


< sup 

0 <e<E O ,raeS JV - 1 


k{n)L ■ n 

c*(n) 


< Too. 


By the parabolic Harnack inequality for strong solutions (see [18, Chapter VII] for instance), we 
get 

V(t 0 ,Xo) £RxR N , max v(t,x) < C2V (t 0 + ^^-r^,xo + k(n)L] , (31) 

*o —1<*<*0, |a=—mo | <1 \ C*(n) ) 


for some C 2 > 0 which is also independent on to, xq, small £ > 0 and n £ S N 1 . 

Combining (30), (31) and the space-time periodicity of the traveling wave, we get 

V(i 0 ,x 0 ) £ R x R n , \d t v(t 0 ,x 0 )\ + \^ x v(t 0 , x 0 )| < C 3 v{t 0 ,x 0 ), 
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with C 3 = C\C 2 ■ Now recall that U*(z,x-,n) = u* ^ x c ^ n [ , x'j ■ Thus 


\d zz u;\ = -^\d t v\ < ^ d z u 


c*{n) 2 


c*(n) 2 c*(n) 


\V x d z U:\ < 


-1 


2 *{n)'- 


dtv n — 


£(«) 


\7 x v 


< 


c*(n) 2 c*(n) 


C 3 v = - 


2 (") 


+ 1 C 3 d z U* 


Since n = info< E < eo inf ne gjv-i c*(n) > 0, this proves the lemma. 


□ 


Let us now go back to the proof of Theorem 4.1. Proceed by contradiction and assume that 
Coo > c*(n) + <5 for some <5 > 0 . We plug U* n k ) into equation (7) satisfied by £/*(•,-;n) and, 
thanks to the above lemma, get 


> 

> 


di v x (AV x U* k ) + (n ■ An) d zz U* k + div x (And z U* k ) +d z (n- AX7 X U*J 
+ q ■ S7 x U* k +(q-n ) d z U* ek + c*{n)d z U* Ek + f(x, U* h ) 

(n ■ An — nk ■ An k )d zz U* k + div x ((An - An k ) x d z U* k ) 

+d z ((n - n k ) ■ A\7 x U* k ) + {q ■ (n - n k ))d z U* k 
+{c{n) - c* Sk (n k ))d z U: k +f(x,U: k ) - f eh {x,U* ek ) 


Aa 2 C\n - n k \ + div x A{n - n k ) + |g||n - n k \ 


5_ 

2 





> 0 , 


(32) 


provided k is large enough, and where a 2 > 0 comes from (3). We now use the sliding method. 
From the asymptotics 

U* k (+oo,--,n k ) = -£k < 0 = C/*(+oo, •; n), 

U* k (-oo,-;n k ) = 1 - y < 1 = U*(+oo, •;n), 

one can define 


t 0 := inf{r : U* k (z + r, x\ n k ) < U*(z, x; n),\/z G R, Va; G 1 N } € R. 

Then, using again the asymptotics as well as the periodicity with respect to x of any pulsating 
wave, there is some first touching point (zo,xq) G R x R w such that 

U e k ( z o + T 0 ,x 0 -,n k ) = U*{z 0 ,x 0 -,n), and U* k (• + r 0 , •; n k ) < U*(-, ■; n). 

Substracting the equation (7) satisfied by U*(z,x-,n) to the inequality (32) satisfied by U* (z + 
To,x;n k ) above, and estimating it at point (^Oj^o); we get that 

£ 

0 > -~d z U* k {z 0 + T 0 ,x 0 ;n k ) > 0, 

a contradiction. Hence, Coo = c*(n), and the convergence of c*(n) to c*(n) is uniform. □ 

Remark 4.3 (On the convergence of profiles). The argument above also shows that the ignition 
traveling waves converge locally uniformly, up to a subsequence, to a traveling wave with minimal 
speed of the monostable equation. Proceeding as in Section 3 and thanks to the monotonicity of 
traveling waves, one can check that this convergence is actually uniform in time and space. In 
particular, they do not flatten as the parameter e —> 0. However, as the uniqueness of the monos¬ 
table traveling wave with minimal speed is not known [If], we cannot conclude on the convergence 
of the whole sequence. 
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5 The uniform lower spreading 

In this section and the next, we will prove Theorem 2.7 under the monostable assumption. The 
easier ignition case will be dealt with in the last section. 

We begin here with the uniform lower spreading property (13) of Theorem 2.7. The argument 
again relies on the approximation from below by an ignition type problem, and follow the footsteps 
of the proof of Theorem 4.1. 

Proof of (13). Recall that f e (x 7 u) is an ignition type nonlinearity which approximates f(x 7 u) 
from below as e —> 0. We still denote u*{t,x) = Uf{x ■ n — c*{n)t 7 x\n) the unique ignition 
pulsating traveling wave of (28) in the direction n, normalized by min xeR w f/*(0,x;n) = 

As / e < /, it is clear that u* E is a subsolution of (1), whose speed is arbitrary close to c*{n) 
as £ —> 0 thanks to Theorem 4.1. This leads back to Weinberger’s result [24], namely the fact 
that for any planar-like initial datum in the n direction, the solution of (1) spreads with speed 
“at least” c*(n) in the n direction. 

Let us now make this spreading property uniform with respect to the family of solutions 
(un) n es N ~ 1 > as stated in Theorem 2.7. In the following /x and K are as in assumption (11) 
(monostable case). Let a > 0 and S > 0 be given. In view of assumption (12) and the comparison 
principle we have u n {t,x) < 1. Hence to prove (13), we need to find r > 0 so that 


inf inf u n (t,x) >1 — 5, 

1 x-n<(c*(n)—ai)t 


(33) 


holds for all t > t. 

In view of Theorem 4.1, we can fix e > 0 small enough so that, for all n £ 

c*(n)>c*(n)-|. (34) 

We then claim that one can find some t e > 0 such that 

u n {t e , x) > 1 — |, (35) 

for all n £ S^ -1 and all x such that x ■ n < —K. We insist on the fact that t e does not depend on 
n £ E N ~ X . To prove (35), let us define 

S = {x £ R jV : x ■ n < c*{n) for all n £ S^ -1 }. 

We know from Theorem 2.5 that the mapping n <— > c*(n) is positive and continuous, hence S has 
nonempty interior. It is then known (see Theorem 2.3 in [24], as well as Remark 2.4 above) that 
for compactly supported initial data “with large enough support”, the associated solution of (1) 
converges locally uniformly to 1 as t — > +oo (in fact, even uniformly on the expanding sets tS' 
for any subset S' of the interior of S; also, under the additional assumption that 0 is linearly 
unstable with respect to the periodic problem, this is even true for any non trivial and compactly 
supported initial datum, regardless of its size [7], [6]). More precisely, let u_ r be the solution of 
(1) associated with the initial datum uo t n{x) = /x x Xb r (x), where R is a large but fixed positive 
constant (depending on /x) which we can assume to be larger than 2 \/N max^ Li. Here Br denotes 
the ball of radius R centered at the origin. Then ur converges locally uniformly to 1 as t —> +00. 
In particular, 

£ 

UR{t e ,x) > 1 - -, 

for some t E > 0 and all x £ B 2 r. Besides, for Xq £ H?L 1 Li'Z such that xq ■ n < —K — R , we have 
— thanks to (11)— that u n (0 7 x + £0) > ur{ 0,x). Then, by the comparison principle, 

£ 

Vx £ B 2 r, u n (t E ,x + x 0 ) > u R (t e ,x) > 1 - 
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Since R > 2 y/N max^ Li, for all x • n < —K, there exists xq £ n^ =1 LjZ such that xq n < —K — R 
and x £ B2r(xo). Thus, we obtain u n {t E ,x) > 1 — §, for all n £ and x ■ n < —K, that is 

claim (35). 

Now, recall that U*(-,-;n) is the pulsating traveling wave of equation (28) in the direction n, 
connecting 1 — | to —e. Hence, it follows from (35) that, for any n £ S^' 1 , one can find some 
shift Z n such that 

Unite, x ) > U*(x ■ n— c*{n)t e + Z n ,x;n). (36) 

Actually, it suffices to select 

Z n := min{z £ R : min U*(—K — c*in)t e + z, x; n) < 0} £ (0, oo). 

Moreover, from the uniform continuity of ignition traveling waves with respect to the direction, 
namely Theorem 2.6, it is straightforward that the family {U*{z, x; n)) ng §jv-i converges to — e as 
z — > +oo uniformly with respect to n £ E> N_1 . Therefore, we can also define the bounded real 
number Z := sup n6 §iv-i Z n £ (0, oo), so that (36) is improved to 

Vn £ S^ -1 , u n {t E , x) > U*{x ■ n— c* e {n)t E + Z, x; n). 

Then we can apply the parabolic comparison principle to get 

Vf > f E ,Vx £ R N ,Vn £ S^ -1 , u n {t,x) > U*{x ■ n - c* e {n)t + Z,x;n). (37) 

Therefore it follows from (34), (37) and the monotonicity of the front that 

u n {t,x) > U* (^-^t + Z,x-nj , (38) 

for all n £ all t > t E and all x such that x ■ n < (c*(n) — a)t. Using again the uniform 

continuity of ignition traveling waves with respect to the direction, namely Theorem 2.6, one can 
find some shift Z' > 0 such that, for all n £ S N_1 , 

z< -Z' =$■ U*{z,x-,n) > 1-e. (39) 

Up to decreasing e, we can assume that e < S without loss of generality. Now choose r >t E such 
that — j t + Z < Z'. Then, we get from (38) and (39) that 

u n {t,x) > 1 - S, 

for all n £ S N ~ 1 , t>r and x such that x ■ n < (c*(n) — a)t. We have thus proved (33), and hence 
(13). □ 

6 The uniform upper spreading 

We conclude here the proof of Theorem 2.7 (monostable case), by proving the uniform upper 
spreading (14) in subsection 6.1. Then in subsection 6.2 we again prove (14) — together with the 
uniform lower spreading property (13)— when assumption (12) is relaxed. 

6.1 Proof of (14) 

We begin by proving some kind of uniform steepness of the monostable minimal waves, which in 
turn will easily imply (14). 

Proposition 6.1 (Steepness of critical waves). Assume that f is of the spatially periodic monos¬ 
table type, i.e. f satisfies (2) and Assumption 1.1. 

Let u*(t,x;n) = U*{x ■ n — c* (n)t, x;n ) be a family of increasing in time pulsating traveling 
waves of (1), with minimal speed c*(n ) in each direction n £ S^ -1 , normalized by U*{ 0,0;n) = 
Then, the asymptotics U*(—oo,x;n) = 1, U*{oo,x\n ) = 0 are uniform with respect to n £ 
gjv- 1 . Moreover, for any K > 0, we have inf ng §.w-i infp|< A ' inf^ggN — d z U*{z, x\ n) > 0 and 
inf/jgs^- 1 inf|z|<ic inf^gRjv U*{z,x\n) > 0. 
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Remark 6.2 (Lack of uniqueness). Such a family of traveling waves is always known to exist. 
However, the uniqueness of the traveling wave with minimal speed in each direction is not known. 
We shall prove that any sequence of increasing in time traveling waves with minimal speed in the 
directions Uk —> n converges, up to extraction of a subsequence, to an increasing in time traveling 
wave with minimal speed in the direction n, as we did in the ignition case. The proposition then 
easily follows, but the lack of uniqueness is the reason we state this result in a slightly different 
way. 

Proof. Proceeding as explained in the above remark, choose some sequence nk — > n £ As 

before, one can extract a subsequence such that u*(-, •; rik) converges locally uniformly to a solution 
Uoo of (1). By the continuity of the speeds c*(n) with respect to n, as proved in Theorem 2.5, the 
function u 0 0 also satisfies 

l ■ n 
c* ( n) ’ 

Moreover, it is nondecreasing in time, hence increasing in time by applying the strong maximum 
principle to dtUoa. In particular, it converges to two spatially periodic stationary solutions as 
t —> ±oo which, as before and thanks to the monostable assumption, must be 0 and 1. As 
announced, u 0 0 is an increasing in time traveling wave with minimal speed in the direction n. 
Reasoning by contradiction, it is now straightforward to prove Proposition 6.1. □ 

Proof of (14). First, from Proposition 6.1 above, and hypotheses (10)—(12), one can find some 
shift K i > 0 large enough so that, for any n £ S^" 1 , uo, n (x) < U*(x -n — Ki,x\ n ). Thus, by the 
comparison principle, 

u n (t, x) < U*(x ■ n — c*(n)t — K\,x; n). 

For any a > 0 and S > 0, let r be such that U*(ar — K\,x\ n) < 5 , for all n £ and x £ R^, 
which is again made possible by Proposition 6.1. Then (14) immediately follows. □ 

6.2 Relaxing assumption (12) 

We here consider the case when the family (wo,n)neS iV - 1 does not necessarily satisfy (12), but is 
only uniformly bounded: there is M > 0 such that 

Vx£R N ,\/n£§ N ~ 1 , u 0 , n (x)<M. (40) 

We prove that, in such a situation, the uniform lower and upper spreading properties (13) and 
(14) remain true if we make the following additional assumptions on the behavior of /, and in 
particular on its behavior above the stationary state p. 

Assumption 6.3 (Additional assumptions). (i) There is <f>(t,x ) a solution of (1) such that 

<j)(0 , •) > M, and <j>(t,x) converges uniformly to 1 as t -± ±00. 

(ii) The steady state 0 of (1) is linearly unstable with respect to periodic perturbations. 

(Hi) There exists some p > 0 such that f(x,u) is nonincreasing with respect to u in the set 
R N x (l-p,l±p). 

The first part of this assumption holds true, for instance, if f(x, s) < 0 for all x £ R^ and 
s > 1. As we will see below, the second part can be expressed in terms of some principal eigenvalue 
problem, and holds true as soon as d u f(x, 0) is positive on a non empty set. The last part is a 
natural extension of (Hi) of Assumption 1.1. 

Combining (33), whose proof does not require assumption (12), and a comparison of the 
solutions (u n ) ne §N-i with (f> given by the above assumption, it is clear that the lower spreading 
property (13) still holds true. In the sequel, we prove the upper spreading property (14). We start 
with the following proposition, whose proof is identical to that of Proposition 6.1 and does not 
require Assumption 6.3. 



VI u 00 (t,x)=u 00 it + 


i =1 
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Proposition 6.4 (Steepness of noncritical waves). Assume that f is of the spatially periodic 


monostable type, i.e. f satisfies (2) and Assumption 1.1. 

For any a > 0, let u a (t , x; n) = U a (x ■ n — (c*(n) + a)t, x; n) be a family of increasing in time 
pulsating traveling waves of ( 1 ), in direction n, with speed c*(n)+a, normalized by U a ( 0, 0; n) = . 

Then, the asymptotics U a (—oo,x;n) = 1, U a (oo,x\n) = 0 are uniform with respect to n £ 
gw-i. Moreover, for any K > 0, we have inf„ 6 §jv-i infu|<^ inf x£R jv — d z U a (z, x; n) > 0 and 
inf„ e s N -i inf|z|<x inf^gRJv U a (z,x;n ) > 0 

We now turn to the proof of the upper spreading property (14), which relies on the construction 
of a suitable family of supersolutions that were already used in [14] (following an idea of [11]). 

Proof of (14). Let a > 0 and S > 0 be given. We need to find r > 0 so that estimate (14) holds 
for all t > t. 

First, we need to introduce some notations, and some well known results (see [25], [4], [14] 
among others). We begin with the principal eigenvalue problem 


-L 0t n,\(t>n,\ = Mo (n, A )(/>„, A in 1^, 

4>n, a is periodic, <j> n> \ > 0, U^aIIoo = 1, 


(41) 


where 


Lo, n> \<j) = div (AVcf) + A 2 (n • An)<t> — A(div(An^>) + n ■ A\/f>) + q ■ Vfi — A (q ■ n)f> + d u f(x, 0 )<j>. 
This arises, similarly as in Section 3, when looking for moving exponential solutions of the type 


e -A 

\{x) of the linearized problem around 0. Such solutions exist if and only if 



which is well-defined thanks to the linear instability of 0 , which reads as /xq = p,o(n,0) < 0 . 


Moreover, it is known that c*(n) > c* in {ri) [13]. We introduce A(n) the smallest positive solution 
of —po(n, A) = (c*(n) + A. It is standard that fio(n, A) is continuous with respect to n and, as 
it is known to be concave, A(n) is also continuous with respect to n. In particular 


0 < min A(n) < max A(n) < +oo. 


n£neS ^" 1 


Let some smooth and nonincreasing x( z ) be such that 



and define, for s > 0 (a shift to be fixed later), 


4>(f, x) = $ s (t, x; n) := x(^ s ) + (1 - x{£.s))<t>n,\{n) (x)e x{n)t 


where 



Note that $ is nonnegative and, along with its derivatives, is bounded uniformly with respect to 
n and s. 

Let us now define various positive constants. Choose 0 < p < 6 small enough so that 



(42) 


( 43 ) 
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Now, by Proposition 6.4, there is K > 1 large enough such that, for all n £ 8^ 1 , x £ 1^, 


£ > K => 0 < (£, x; n) < -, 


£<-K =>!--< U<±(£,x\n) < 1. 


( 44 ) 


Then, by Proposition 6.4 again, we have 


7 := inf inf 

ries "- 1 |z|<jf,aeR JV 


-d z U°(z,x-,n) > 0. 


(45) 


Last, we define 

rj a'y 

61 := 2p|U ’ 62 := 4(||9t$||oo + ||div(^4V4 > )|| 00 + || 9 • V$||oo + ||$||oo||a u /|| i =c (RWx(0il+3)) ) 

and 

e := min (ei, e 2 ) > 0. (46) 


Now, we are going to show that 

v(t, x) = v s (t, x;n) := [x ■ n — (c*(n) + —^ t — s, x\ nj + e<f >(t, x) = U » (£ s , x\ n) + e<f>(t, x) 

is a supersolution of the monostable equation (1). Straightforward computations and the mean 
value Theorem yield 


£[x](f, x) := d t v(t,x) — div(A(x)Vu(f, x)) — q(x) ■ X7v(t,x) — f(x,v(t,x)) 

= e [d t $(t, x) - div(A(x)V4>(f, x)) - q(x) ■ V<f>(i, x) - $(t, x)d u f(x, 0{t, a;))] 
-^dzU%{t; s ,x-,n), 

for some 

£/*(£ s ,x;n) < 8(t,x) < U°(£ a ,x;n) + e®(t,x). 

We distinguish three regions, depending on the values of £ s . 

First, if |£ s | < K, the nonnegativity of C[v](t,x) is obtained thanks to —jd z U<±(t; s ,x-,n) > fj 
by (45) and the definition of e in (46). 

Next, if £ s > K , then 4>(t, x) reduces to <j> n ,\(n){x)e~ x ^ n ^ a and, dropping —jd z U<*((; s ,x-,n) 
which is positive, we arrive at 


-C[v](t,x) > 

> 


A(n) (c*(n) + ^) + Mo (n, A (n)) + d u f(x, 0) - d u f(x, 0(t , x)) 
(-j-A(n) + d u f(x, 0) - d u f(x, 0{t, x))) 0 n>A (n) (*> 


(t>n,\{n)(x)e A(n)5s 


o-A («)£» 


But, when £ s > K , (44) and e < ei imply 0 < 8(t,x) < rj, and the nonnegativity of C[v\(t,x) is 
obtained thanks to (42). 

Last, we consider the case where £ s < —K, so that $(t, x) reduces to 1. Hence 


^C[v\(t,x) > —d u f(x, 0(t, x)). 

But, when < —A', (44) and e < ei imply 1 — 77 < 0(t,x) < 1 + 77, and the nonnegativity of 
C[v](t,x) is obtained thanks to (43). Hence, v s (t,x;n) is a supersolution of (1). 

Thanks to (40), we get by the comparison principle that, for all n £ S^ -1 , all t > 0, all 
x £ R-^, u n (t,x) < (j>(t,x), where (f> is given by Assumption 6.3. Now choose T > 0 such that 
(j>(T, x) < 1 + |, and get that 


Vn £ 1 , Vx £ R^, u n (T,x) < 1 + ^. 


(47) 
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Using the comparison principle and a computation identical to that of (17), we get that, for 
any large A > 0, there is C > 0 - independent on n thanks to (10) and (40)— such that 

Vn € S N ~\ Vi > 0, Vx G R n , u n (t,x) < Ce" A(a: ' n - 2oiAt) . 

In particular u n (T , •) decays faster than any exponential as x • n —> +oo, namely 

VA > 0, u n (T, x)e Xx ' n —> 0 as x ■ n —> +oo, uniformly w.r.t. n G S^ -1 . (48) 


Observe that, for all s > 0, 

Vn G § N ~ 1 , Vx-n > ^c*(n) + ^ T+s+1, v s {T,x\n) > e<f) n> \( n )( x )e~ X( ' n ' >x ' n > e ye _Am “ xX '", (49) 


where 7 := min„ g §jv-i min xgH iv (j) n t \( n ){x) > 0 and X ma x '■= max„ 6 gw-i A(n) < 00 (recall that 
n A(n) is continuous and so is (n, A) Now, select A > 1 large enough so that, for all 

s > 0 , 

Vn G Vx ■ n < (c*(n) + ^ j T + s — A, v s (T, x; n) > 1 + (50) 

which is possible thanks to Proposition 6.4, and more precisely the uniform with respect to n 
asymptotics of U<*(z,x-,n) as 0 —> — 00 . Proposition 6.4 also enables to define 


k := 


inf inf inf 

raes N ~ 1 -A<z<1i£K» 


U% (z,x;n) > 0, 


so that, for all s > 0, 

Vn G § N ~ 1 , V ^c*(n) + —^ T + s — A < x ■ n < ^c*(n) + — ^ T + s + 1, v s (T , x; n) > k. (51) 
In view of (48), we can now select a large enough shift sq > A so that 

Vn G Vx ■ n > (c*{n) + ^ T + so — A, u n (T,x ) < minjey, re}e _AmaxX ' rl . 

Combining this with (47), (49), (50), (51), we have that, for all n G and x G R^, 

u n (T,x) < v So (T, 2; n). 


Then, by the comparison principle, for all t > T, x G R ff , n G § w 1 , 

0 < u n {t, x) < U « (x ■ n — (c*{n) + — j t — so, x\ nj + e4>(i, x). 
Hence, when x ■ n > (c*(n) + a)t, we have, since e^Hoo < f < |, that 

0 < u n (t , x) <U« (^t - s 0 , X] nj + | < 5, 


as soon as t > r, where r > 0 is large enough (again independently on n by Proposition 6.4). This 
proves (14). □ 


7 The uniform spreading: the ignition case 

For the sake of completeness, we give here the main steps to prove Theorem 2.7 in the (simpler) 
ignition case. We will see that it follows from the continuity of ignition waves, Theorem 2.6, 
together with the standard idea explained in Remark 2.4. We will briefly sketch at the end of this 
section how the hypothesis (12) can again be relaxed. 
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Proof of Theorem 2.1 in the ignition case. First, the proof of the uniform upper spreading (14) 
is the same as that of subsection 6.1 in the monostable case, using Theorem 2.6 (continuity of 
ignition waves) instead of Proposition 6.1 (steepness of critical waves). 

Let us now prove the uniform lower spreading (13). Let a > 0 and <5 > 0 be given. We 
may reduce S without loss of generality, and assume that S < p where p is given by part (Hi) of 
Assumption 1.2. Using the same arguments as in the proof of (35), we get the existence of some 
time ts > 0 such that 

£ 

U n (ts, X) > 1 - 

for all n £ and all x such that x ■ n < —K. Now let, as usual, U*(x ■ n — c*(n)t,x\n) be 

the unique ignition pulsating wave in the direction n, normalized by min ieH jv U*(0,x;n) = 
Thanks to the inequality above and the continuity of the mapping n4 [/*(•,•; n) with respect to 
the uniform topology, it is clear that there exists some shift Z > 0 such that, for all n £ SA -1 , 

U*(x ■ n + Z,x;n) — — < u n (ts,x), \/x £ R N . 

We then check that u(t, x) := U*{x ■ n + Z — (c*(n) — f)t,x-,n ) — § is a subsolution of (1). 
Indeed, 

d t u — div (A(x)Vu) — q{x) ■ Vu — f{x,u) = ^ d z U* + f(x, U*) — f{x,u). 

Assume first that u < 9 — f. Then f{x,u ) = f(x,U*(x ■ n + Z — ( c*(n ) — f)t,x;n ) = 0, which 
together with the monotonicity of U* with respect to its first variable, gives the wanted inequality. 
Assume then that u > 1 — p. Then, by the monotonicity of / with respect to u in the range 
[1 — p, 1], we again obtain the wanted inequality. 

It remains to prove that u is a subsolution when 9 — | <u< 1 — p or, equivalently, when 
6 < U*(x ■ n — (c*(n) — j)t + Z,x\n) < 1 — p + |. Recall first that 1 — p+| < 1 — §■ Using 
again the continuity of the ignition wave with respect to the direction, we have that there exists 
some R > 0 such that, for all n £ $ N ~ 1 , 

z>R => U*(z + Z,x',n) < 9, z <—R => U*(z + Z, x;n) > 1 — 

and, furthermore, 

max max max d z U*(z + Z,x;n) < 0. 

neS N ~ 1 |z|<fl cceK^ 

Up to reducing <5 again, we may assume that 

max max max d z U*(z + Z,x'n) < -, 

nes"- 1 |«|<r xeR N a 

where M is a Lipschitz constant of /. Therefore, when 9— | <u< 1 — p, then \x-n— (c*(n)— ^)t\ < 
R and 


dtu — div (A(x) Vu) — q(x) ■ Vu — f(x, u ) 

= f- d z U* + f(x,U*)-f(x,u )I 

< ^d z U* + M 5 -< 0, 

that is the wanted inequality. 

We can therefore apply the comparison principle and conclude that 

U* (x ■ n + Z - ( c*(n ) - C -^)t,x;nj - | < u n (t s +t,x), 
for all x £ and t > 0. 
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Noting that there exists some other shift Z' > 0 such that, for all n £ 1 , 

z < —Z' => U* (z, x\ n) — — > 1 — 5, 

we get the uniform lower spreading (13) as in the end of Section 5. □ 

Relaxing hypothesis (12). In order to relax (12), assume now that / satisfies Assumption 1.2 and 
parts (*) and (Hi) of Assumption 6.3. As above, one can then show that U*(x ■ n — ( c*(n) + 
Tj)t, x\ n) + | is a supersolution of (1). Then, as in Section 6.2, one can find some time T and some 
shift s 0 such that, for all n, the solution u n (T , x) lies below U* (x ■ n + sq — (c* (n) + f)T, x; n) +1 in 
the whole space. It is then straightforward to obtain the wanted uniform upper spreading (14). □ 
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